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ABSTRACT 

We report on the discovery of three bright, strongly-lensed objects behind 
Abell 1703 and CL0024+16 from a dropout search over 25 square arcminutes of 
deep NICMOS data, with deep ACS optical coverage. They are undetected in the 
deep ACS images below 8500 A and have clear detections in the J and H bands. 
Fits to the ACS, NICMOS and IRAC data yield robust photometric redshifts in 
the range z ~ 6 — 7 and largely rule out the possibility that they are low-redshift 
interlopers. All three objects are extended, and resolved into a pair of bright 
knots. The bright i-band dropout in Abell 1703 has an if-band AB magnitude 
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of 23.9, which makes it one of the brightest known galaxy candidates at z > 5.5. 
Our model fits suggest a young, massive galaxy only ~ 60 million years old 
with a mass of - W^Mq. The dropout galaxy candidates behind CL0024+16 
are separated by 2.5" (~ 2 kpc in the source plane), and have if-band AB 
magnitudes of 25.0 and 25.6. Lensing models of CL0024+16 suggest that the 
objects have comparable intrinsic magnitudes of AB ~ 27.3, approximately one 
magnitude fainter than L* at 2; ~ 6.5. Their similar redshifts, spectral energy 
distribution, and luminosities, coupled with their very close proximity on the sky, 
suggest that they are spatially associated, and plausibly are physically bound. 
Combining this sample with two previously reported, similarly magnified galaxy 
candidates at z ~ 6 — 8, we find that complex systems with dual nuclei may be 
a common feature of high-redshift galaxies. 

Subject headings: cluster: general - cosmology: observation - galaxies: high- 
redshift 



INTRODUCTION 



Deep observations and surveys with the Hubble Space Telescope (HST), suc h as GOODS 



( lOesch et al. 



Giavalisco et al.ll2004J : [Dickinson et al.ll2004f ) , the HUDF flBeckwith et a 



20071 ). and the HUDF parallel observations (IBouwens et al. 



2006), the HUDF05 



2004d) have greatly 



enriched our knowledge of the distant Universe, enabling the discovery of more than 600 
^775— band dropouts (i-dropout; 



2;-dropout; iBouwens et al. 



Bouwens et al.ll2006l . 120071 ). and eight ^sso— band dropouts 
20081 ) which are believed to be galaxies at 2; ~ 6 and 7, respec- 
tively. An L* galaxy of M = —20.3 at 2; ~ 6 exhibits an apparent magnitude ~ 26.8 (AB, 
and hereafter); consequently, the majority of 2; > 5.5 galaxies are fainter than AB=26. 

Bright high-redshift galaxies are rare. Ten i-dropouts brighter than AB = 25.5 were 
found in ~ 400 square arcminutes of GOODS, HUDF, and HUDF parallel fields ( IBouwens et al, 



20061. 



200 



i 



( jBunker et aL 



wo of them have been con firmed to be high-redshift objects spectroscopically 



2003 



Stanwav et al. 



arcminutes (IShimasaku et al 



2005 



2004r). In a survey of th e Subaru Deep Field of 767 square 



Taniguchi et al.ll2005l ). 18 i-dropout galaxies were spec- 
troscopically confirmed to be at z ~ 5.7, and nine at 2 ~ 6.6, but the brightest among them 
is only at magnitude of 25.4. In a recent s urvey that covers o ne square degree in the field 
of the Subaru/XMM-Newton Deep Survey (jOuchi et al.ll2008a ). 17 Lya emitters (LAEs) at 
z ~ 5.7 were found, of which only one is brighter than AB=25.5 in the z band. At least one 
bright LAE has been found at z ~ 6.6 (z-band magnitude 25.4, narrow band (921nm) mag- 
nitude 23.5; lOuchi et al.ll2008al ). Combining the results of these three surveys, the brightest 
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z ~ 6 candidate over 5000 square arcminutes is at AB=24.5 (in broad bands - LAEs are 
intrinsically brighter in a narrow band). 

Gravitational lensing by massive galaxy clusters enables the detection of high-redshift 
galaxies that are fainter than can be found even in very deep surveys. In addition, lensed 



(Hu et al. 


2002; 


Kneib et al. 


2004 


; 1 


„, 

^ello et al. 


2004 




Schaerer & Pelld 


2005 




Richard et al. 


2006; 


Frve et al. 


2008; 


Bradlev et al. 


20081; 


Bouwens et al. 


2009al; 


Richard et al.l 


20081). The 



number of z > 6 candidates reported, however, has varied considerably depending on how 
aggressive the selection criteria are (and conse quently how reliable the detected sources are: 
see, e.g., discussions in iBouwens et al.l l2009d ). The HST images of clusters are the best 
in terms of spatial resolution and photometric accuracy, thus permitting the most reliable 
selections of high-re dshift galaxie s . On e salient example is the three i-dropout arclets in 
Abell 2218 found by iKneib et al.l (120041 ). These arclets represe n t diffe rent views (multiple 
images) of the same source at z = 6.75 ± 0.1. iBouwens et al.l (120031 ) found an i-dropout 
object in the field of cluster RDS1252— 29 as bright as AB=24.2, and follow-up spectroscopy 
( iDow-Hygelund et al.ll2005l ) confirmed it as a Lyma n-break galaxy at z = 5.515 and revealed 
stellar absorption lines in the continuum. Recently, iBradley et al.l (120081 ) discovered a galaxy 
candidate at z ^ 7.6 in the field of Abell 1689 that is spatially resolved and has an observed 
if-band magnitude of 24.7. The cluster lensing magnifies this galaxy (Al689-zDl) by a 
factor of ~ 10. 

Taking advantage of th e NICM OS data of 21 arcmin^ over 11 galaxy clusters (including 
Abell 1689) . IBouwens et al.l (l2009al ) searched for strongly lensed 2;-dropout galaxies at z > 7 
and found that Al689-zDl is the only robust candidate in these fields. This suggests that 
bright high-redshift galaxies may be rare in cluster fields, even given the smaller field size of 
cluster regions. Since the amount of NICMOS data available over cluster fields is still very 
limited, there is reason to look at the surface density of such objects over larger areas and 
at somewhat lower redshifts, in order to overcome the effect of small number statistics. 

The purpose of this paper is to report on the results of a search for bright (if-band 
magnitu de < 25.6), high l y mag nified i- and z-dropout galaxies over a larger area than was 
used by IBouwens et al.l ( l2009a[ ). Such sources provide us with a unique opportunity to 
examine the properties of lower-luminosity galaxies in detail. The lower- luminosity galaxy 
population is of great interest because it appears to dominate th e UV luminosity density , 



Yoshida et a. 



SFR density, and production of ionizing 



Reddy et al. 



20061; 



20081; 



Beckwith et al 



Frye et al 



20081; 



( An < 912A) photons 



2006 



Bouwens et al. 2006 



Chen et al. 



Yan fc Windhorst 



2007 



Oesch et al. 



2004 



2007 



2009! 



]). In particular, we consider ~ 4 square 



arcminutes of additional search area over CL0024-(-16 and CL1056— 03 obtained in a recent 
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NICMOS program HST-GO10874 (PI W. Zheng). These areas were not included in the 
z-dropout search of iBouwens et ahl (l2009d ). 

The plan for this paper is as follows. In §2, we describe the observational data. In §3, 
we outline our search criteria and summarize the results of our search for bright z ^ 6 — 7 
galaxies. In §4, we use stellar-population modeling and the available lensing models to 
remark on the physical properties of the high-redshift galaxy candidates we are examining, 
and in S5, we summarize our results. 



2. OBSERVATIONS 

A summary of the HST ACS, WFPC2, and NICMOS imaging observations used here 
to search for 2 ~ 6 — 7 i- or 2;-dropout galaxies is given in Table [H Of particular interest 
are new NICMOS observations acquired over CL0024+16 (two NICMOS NIC3 pointings) 
and CL1056-03 (three NICMOS NIC3 pointings) from the program HST-GO10874. These 
observations were taken to better quantify the volume density of faint galaxies at z ~ 7. 
Each NICMOS image is only 52" x 52" in size and covers just a portion of the total area 
with deep ACS data over these two clusters. 

The NICM0S/NIC3 images of CL0024+16 were obtained in 2007 July, with exposures 
of 9.4 ksec in both the FllOW and F160W band. The NICM0S/NIC3 images of CL1056-03 
were obtained in 2007 January and December, and 2008 January, with exposures of 6.6-8 ksec 
in both the FllOW and F160W band. We used a custom pipeline written at the Unive rsity 



of California Santa Cruz to process the NICMOS images ("nicred"; iMagee et al.l 120071 ) into 



a final mosaic. The 5a limiting magnitudes for these data are approximately 26.7. 

CL0024+16 was observed in 2004 November in six broad optical bands: F435W, F475W, 
F555W, F625W, F775W and F850LP, including 16.3 ksec in the z-band and 10.1 ksec in 
the i-band. CL1056-03 was observed in 2004 March in the F606W, F775W and F850LP 
bands, including 17.7 ksec in the 2;-band and 17.3 ksec in the z- band. All the ACS sc ience 



images over these fields were processed with our pipeline APSIS ( jBlakeslee et al.ll2003l ) and 
registered with the available NICMOS data. 

We retrieved archival Spitzer/IRAC images of the field of Abell 1703 and CL0024-I-16, 
observed in 2007 January and 2003 December, respectively. The mosaic images are in the 
form of Post-BCD (Basic Calibrated Data), produced by the Spitzer pipeline calibration 
MOPEX. For Abell 1703, the exposure times were 6.2 ksec for both the 3.6 and 4.5 fim 
bands, reaching a 5a limiting magnitude of ~ 24.7. For CL0024+16, the exposure times 
were 2.7 ksec in the 3.6yum band and 4 ksec in the 4.5 /xm band, reaching a 5a limiting 
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magnitude of ~ 24.5. 



We also searched for bright 2; ~ 6 z-dropouts in the cluster data considered by lBouwens et al. 



( l2009d ). who had used these data to search for z- and J-band dropouts at z > 7 and in- 
clude NICMOS observations over 11 high-redshift clusters (Abell 1689, Abell 1703, Abell 
1835, Abell 2218, Abell 2219, Abell 2390, Abell 2667, AC114, CL0024+16, 1E0657-56, and 
MS1358+61). Most of these clusters have deep ACS and WFPC2 optical data to ensure that 
the candidates we select at z ~ 6 as i-dro pouts are not detected in the optical to very faint 



limits. We refer to iBouwens et al.l (j2009al ) for more details on the reductions and processing 
of those data. Since we will be reporting on a bright 2; ~ 6 candidate behind Abell 1703, we 
will summarize its ACS + NICMOS observations as well. NICM0S/NIC3 images of Abell 
1703 were obtained in 2007 July, with exposures of 2.6 ksec in both the FllOW and F160W 
band. The cluster was observed at three other positions in the FllOW band (see Fig. 1). 

Independently as part of another program to search for high-redshift galaxies, we have 
acquired deep IR images of these cluster fields using a wide variety of different ground-based 
telescopes (Gemini-North/NIRI, Magellan/PANIC, and the VLT/ISAAC). Abell 1703 was 
observed with the Gemini-N/NIRI on 2003 June 11, with a total exposure time of 7.6 ksec 
in the J band, and a 5cr limiting magnitude of 23.3. The Magellan/PANIC observations of 
CL0024+16 were made in 2003 October 3 and 5: the J-band observations took 17.8 ksec, 
and the image reaches a 5o" limit of AB~ 26. The infrared images were processed with 
the IRAF/XDIMSUM package. Images taken on each night and with the same instrument 
configuration were sky-subtracted, shifted by the amount determined using the positions of a 
bright point source in dithered images, and combined. We select dropout galaxies when they 
are not detected beyond 2a level, and the color decrements are more than two magnitudes. 
While these ground-based observations are shallower than the NICMOS d ata, they were used 



to tentatively identify the i-dropout galaxies. In a forthcoming paper (jZheng et al.l l2009l ) 
we will be presenting a few of these z-dropout galaxies from these wider-area ground-based 
data. 



3. SELECTION OF BRIGHT i- AND ;z-BAND DROPOUTS 



Source catalogs were produced using SExtractor (IBertin fc Arnouts 1119961 ). We coad- 
ded the available N ICMOS images for the clusters under study to create the detection image 
dSzalav et al]ll999h . We then measured the fluxes on each of the ACS, WFPC2, and NIC- 
MOS images using the same scalable apertures, after PSF-corr ecting the higher resolu tion 
ACS and WFPC2 images to match the NICMOS images. As in Bouwens et al.l f2009ah . we 



measure colors from smaller scalable apertures (Kron factor of 1.2; iKron Ill980l ). typically 
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~ 0.3". We then correct the fluxes measured in these smaller apertures to total magnitudes 
by comparing the light inside a larger scalable aperture (Kron factor of 2.5) of ~ 0.6", on 
the detection image to that inside a smaller aperture. 

Our color selection is based on the magnitudes measure d in small-scalab le aper tures. 
The selection c riteria for i-dropouts are similar to that of iGiavalisco et al.l (120041 ) and 
Bouwens et al.l (120061 ) . but with some variations: 

(^775 - ^850) > (1-3 + 0.9(2:850 " ^lio)) 

or 

(«775 - 2;850) > 2.0 A (2:850 - ^iio) < 0.8 

where A represents the logical AND symbol. In addition, we require that objects are not de- 
tected (< 2a) in all bluer bands to exclude lower-redshift interlopers. Also, selected dropout 
galaxies must have a star/galaxy index < 0.8 and an image dimension greater than three 
pixels to exclude stars and artifacts. 



The selection criterion for the 2-dropouts (IBouwens et al.ll2008l ) is: 

(^850 - ^iio) > 0.8)A 

(%50 - ^iio) > (0.8 + 0.4(Jno - ^^i6o))A 

(-^110 — Hieo)) < 1.2 

Also, we require that objects are non-detected in all the bands shortward of 8500A. 

In cases of non-detection in a shorter-wavelength band, we define the color decrement by 
setting the flux in that band to be zero, with an error that is equal to the la detection limit. 
The non-detection in the r,g, V bands is important, as some sources are weakly detected in 
some of these bands. We found five other objects that were weakly detected in the shorter- 
wavelength bands. Their large i — z color decrements are probably produced by a Balmer 
decrement instead of a Lya edge, and therefore the objects are at low redshift, e.g. 2 ~ 1. 
Without deep images in the shorter-wavelength bands, these objects may be considered as 
"high- redshift candidates" . This may explain why some other searches yield a considerable 
number of dropout galaxies. 

Finally, we set an if-band magnitude limit of 25.6 to identify only the brightest sources 
in our survey. Applying the above selection criteria across our search fields, we identify three 
new dropout galaxie s in total. One source, called Al703-iDl, is found in Abell 1703 (see also 
Richard et al.ll2009l . who also report on this source) and two other sources, CL0024-iDl and 
CL00 24-zDl, are found over CL0024+16. We also identified the 2 ~ 6 .5 i-dropout in Abe ll 
2218 jKneib et al.l I2OO4I ) and the z ~ 7.6 ^-dropout over Abell 1689 JBradlev et al.ll2008h . 
No bright dropouts are found in CL1056— 03 or any other clusters. Our three new galaxy 
candidates are given in Table [2l along with the previously known sources. The positions of 
these sources in Abell 1703 and CL0024+16 are shown with the red circles in Fig. [T]|2l Their 
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cutout stamps are given in Fig. 

For our high-redshift candidates, we used the available Spitzer/IRAC data to mea- 
sure fluxes at 3.6 and 4.5 fim. Since these data have a considerably larger PSF than the 
ACS/NICMOS data, this poses a challenge for accurate photometry. First, the IRAC PSF 
(~ 1.6" at 3.6/im) is considerably larger than the apertures used in the HST images (typ- 
ically < 0.5" in radius). We used an aperture size of 2.5" diameter, and then applied 
corrections of 0.56 and 0.6 magnitude for the 3.6/im and 4.5yum band, respectively, to ac- 
count for flux outside of the aperture. These corrections were derived for stellar proflles 
(see, e.g., iBradley et al.l l2008l ) . In the fleld of Al703-iDl, there are two nearb y sources of 



comp arable brightness approximately 2" away (Fig. 3). We ran the task Galfit ( jPeng et al. 



20021 ) to deblend the source fluxes with flxed source positions and a standard PSF profile of 



FWHM=1.6". The derived magnitudes are presented in Table 2. 

Al703-iDl is a prominent i-dropout object, with a z-band magnitude AB= 24.2. As 
shown in Fig. 3, it is well detected in the ACS 2;-band and NICMOS J and H bands, but 
not in the i-band and the other bands at shorter wavelengths. The object is detected in 
the Gemini-N image in the J band, with a magnitude of 24.2 ± 0. 4. Our confidence th at 
this source is a galaxy at z ~ 6 rather than a lower-redshift object (IStanway et al.ll2005l ) is 
largely motivated by the source's large color decrement of >2.6, which is twice the threshold 
in a normal selection (1.3 — 1.5). The possibility that this source is a dusty galaxy at lower 
redshift is ruled out as the source is roughly equally bright in the z- and J-bands, and 
entirely undetected in the rgoe- [r) and z-bands, which is clearly a sign of a break in the 
spectral energy distribution (SED) instead of a gradual decrease toward bluer wavelengths. 
The source is resolved into two components (magnitudes 24.3 and 24.9; Fig. 3) separated 
by ~ 0.45". The only other objects with such colors are T-dwarfs, which can be ruled out 
by the fact that the source is extended. The FWHM for each component is ~ 0.4", which 
is four times the ACS resolution. This source is detected at low significance (< 2a) in the 
F435W band (magnitude 27.7 ± 0.5), but not in an overlapping band F475W (> 28.5) and 
F555W and F625W. This weak detection is probably just an upward statistical fiuctuation 
and does not affect the SED fitting results noticeably. 

The two sources in CL0024 pose an interesting question. These sources, CL0024-zDl 
and CL0024-iDl, are separated by only 2.5". As we shall see below they have very similar 
redshifts and SEDs, and this small separation corresponds to just 2 kpc in the source plane. 
Thus they are likely to be spatially associated. This will be discussed further below, but for 
now we discuss their individual properties. 

The source CL0024-iDl is resolved in both the J and H band images. Its NICMOS J- 
band magnitude is 25.1 (Table 2). An independent check comes from our Magellan/PANIC 
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image data which gives 25.8 ± 1.1, which is consistent with the NICMOS measurement. 
The source CL0024-zDl appears to be marginally resolved into two components in the Juq 
band image (Fig. 4), but not in the Hiqq band. In the Magellan images, its magnitude is 
26.3 ± 1.3, which again is consistent with the NICMOS data. The detection in the 4.5/im 
band is marginal, and appears to be consistent with background fluctuations. 



It is interesting to note that Al689-zDl {z ~ 7.6; iBradley et al.ll2008l ) also consists of 
two components separated by ~ 1 kpc in the source plane. The apparent frequency of high- 
redshift galaxies with two bright components is perhaps not surprising, in that we expect 
to see bright star-bursting galaxies imbedded in massive halos that are merging at this very 
early epoch. 

Because of the long time baseline of our ACS, NICMOS, IRAC, and ground-based 
observations, none of these sources are likely to be supernova outbursts or transiting objects. 
In addition, because the sources are detected in multiple bands, it is highly unlikely that 
they are spurious detections. 



4. DISCUSSION 



Current studies of the UV luminosity function (LF) at z ~ 2 — 6 suggest that lower- 
luminosity galaxies provide the dominant contribution to the overall luminosity density in 
the UV, the stellar mass de nsity, the ionizing photon budget necessary f or reionization, 



and t he build-up of metal s (lYan fc Windhorst 



20061 : iBouwens et al.ll2006l . 120071 : iReddv et al. 



2004 



Yoshida et al. 



2008; Chen et al. 



2006; Beckwith et al. 



20091 ). As such, clearly it 



is important to quantify the properties of this faint population. Unfortunately, this can be 
difficult to do in the field, due to the integration times required and due to crowding issues 
in deep IRAC images. This makes the present sample of highly magnified galaxy candidates 
at z ~ 6 — 7 very useful for investigating the properties of lower-luminosity galaxies at high 
redshifts. 



4.1. Intrinsic Luminosity 

To determine the intrinsic properties of the high-redshift galaxy candidates identified 
in our search, it is necessary to model the effect of gravitational lensing on these sources. 
Fortunately, both clusters behind which our galaxy candidates lie have been the subject of 
substantial modeling efforts. Abell 1703 {z = 0.258) is known to have a lar ge number (i.e. 



16) of multiply lensed galaxies and has been the subject of detailed modeling (ILimousin et al. 
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20081). CL0 024+16 (z = 0.396 ) 



the lensing (IKneib et al 



2003 



has also been the subject of several detailed efforts to m odel 



Jee et al.ll2007l : iNataraian et al.ll2007l : IZitrin et al.ll2009f ). We 



adopt the model of IZitrin et al.l (120091 ) as our fiducial lensing model, but compare with 
other models to gauge the approximate uncertainties. As shown in Table 4, we adopt a 
magnification at our source position of 3.1 for Al703-iDl, 6.3 for CL0024-iDl, and 6.2 for 
CL0024-z pi. For CL0024+16, a wide range of magnification exists in th e literatur e . The 
model of iNatarajan et al.l (120071 ) yields a magnification of 2 and that of I Jee et al.l (120071 ) 
predicts a magri i ficatio n of 17-20. The low magnification value probably is arose because 
Natarajan et al.l (120071 ) assume the cluster mass profile as singular isothermal, and they used 
five multiple images of the same s ource near the center (as well as weak lensing constraints). 
The recent work of IZitrin et al.l (120091 ) takes advantage of high-quality, multi-band ACS 
data and identifies 11 sets of multiple images covering a wide range of background redshift 
0.5 < z < 4.5, so that the gradient of mass profile slope can be meaningfully constrained 
and hence more reliable magnification measurements are obtained. While we have adopted 
perhaps the most reliable magnification model, the actual magnification factors in current 
lensing models are sensitive to details in the mass distribution of the foreground cluster and 
thus subject to significant uncertainties. 

Adopting the photometric redshifts for A1703-iDl, CL0024-iDl, and CL0024-zDl that 
we will determine in the next section {i.e., z = 6.0, z = 6.5, and z = 6.6, respectively), we 
estimate that their amplification factors are approximately 3.1, 6.3, and 6.2. respectively. 
This implies intrinsic J-band luminosities for these sources of AB=25.2, 27.1, and 28.0, 
respectively. The two fainter dropout sources in CL0024-I-16 are estimated to have intrinsic 
magnitudes fainter than L* after demagnif ication and would no t likely have been foun d 
in even the NICMOS data to the HUDF JSouwens et al.l 12004^ : [Thompson et allbood l 
Al703-iDl, on the other hand, appears to be quite bright, having an intrinsic magnitude of 
~ 25. 



4.2. Spectral Energy Distribution and Photometric Redshift 



We performed fits to t he multiband photometry of these three sources by using the 



stellar-population models of iBruzual fc Chariot! (12003). For a ll of th e models, we assume 



a simple stellar population (single-burst; SSP) with a ISalpeterl (Il955l ) initial mass function 
with mass cutoffs of 0.1 and 100 Mq. We explored models with both solar {Z = 0.02) 
and subsolar ( Z = .0004 = Zq/50) metallicities. The models use the obscuration law of 



Calzetti et al.l (120001 ) and are correcte d for Lyman-s eries line-blanketing and photoelectric 



absorption following the procedure of iMadaul (119951 ). In the stellar-population model fits. 
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we constrain the stellar age to be less than the age of the universe at the fit redshift [e.g., 
0.8 Gyr at z = 6.5). 

The best-fit stellar-population models are shown in Fig. 013 and the parameters are 
given in Table [31 While the non-detections are shown in these figures as 2 cr upper limits, 
we used the measured fluxes, not the limits, to fit the SEDs. For both solar and subsolar 
metallicities, we found acceptable fits for the models. IRAC photometry places important 
constraints on the stellar-p opulation models as it provides the rest-frame optical magnitudes 



for high- redshift galaxies (iLabbe et al. II2006I ). Using the ratio of UV-optical luminosity 



one can estimate the galaxy age with considerable accuracy. Considering the uncertainties 
associated with the mass models themselves and in potential reddening factors, the real 
errors are likely somewhat higher than those listed in Table 3. 

For A1703-iDl (Fig. 5), the best-fit redshifts are z = 5.9 — 6.0 with an 1 a error 
of 0.1. The best-fit results, with xl = 1-1; sug gest a high-re d shift galaxy with a mass 



of 0.9 — 1.2 X 10^^ Mq and age of 55 — 64 Myr. iRichard et al.l (120091 ) measure a redshift 
oi z = 5.827 from Keck spectroscopy, which is broadly consistent with our fitting results. 
We also ran the fits to photometry with the redshift fixed to 2; = 5.827, and the results 
are presented in Table 3 and Fig. 5. Using the spectroscopic redshift removes one of the 
degrees of freedom, but the modeling results still do show some dependence on the unknown 
metallicity: The fit with solar metallicity produces a satisfactory result, with xl = l-l? but 
that with a low metallicity does not fit as well, with xl = 1-7, suggesting that the metallicity 
may well be high in this very luminous (3 L*) galaxy, or, if the metallicity is low, that the 
low-metallicity models may not match the SED very well in these high-redshift galaxies. 
The low-redshift solutions are inconsistent {xl ~ 10) with the observed photometry. 

For CL0024-iDl (Fig. 6), again we find that the high-redshift solutions {xl = 0.5 — 0.7) 
provide a better fit to the data than the low-redshift solutions {xl = 2.8 — 4.2). Of note, 
the low-redshift SED models are inconsistent with the observed r- and z-band 1 a upper 
limits and the z- and Jno-band photometry. The best-fit redshifts are z = 6.4 — 6.5 with an 
1 a error of 0.1. The stellar masses are 1.5 — 2.1 x lO^M© with ages of 43 — 72 Myr. The 
reddenings derived from the fits to both the solar and subsolar models were Ay = 0.0. 

For CL0024-zDl (Fig. 7), we find the high-redshift solutions {xl ~ 0.1) are only 
modestly better than the low-redshift solutions {xl = 0.6 — 0.9). With the exception of 
the 2;-band photometry, which has large uncertainty, the low-redshift SED models can also 
reproduce the measured fluxes blueward of the break in the SED. We still favor the high- 
redshift solutions {z = 6.5 — 6.6 ± 0.4) over the low-redshift models {z = 1.4 ± 0.6). The 
stellar masses are ~ 2 — 2.5 x lO^M© and the ages ~ 85 Myr. The solar and subsolar models 
have fit reddenings of Ay = 0.0 and 0.2, respectively. 



4.3. Are CL0024-iDl and C10024-zDl Physically Related? 



As noted in §3, the source CL0024-zDl is located only 2.5" from the CL0024-iDl coun- 
terpart. For our fiducial lensing model where the amplification is 6.2 (Table 4), they would 
be separated physically by only 2.2 kpc. The two sources are of similar luminosity, differing 
by only approximately 0.7 magnitude in both the in J- and if-bands. The SED fits suggest 
very similar redshifts, differing only by Az ~ 0.1 , which makes them essentially identical 
within the errors (see Table 3). The formal designation indicates that one is an i-dropout 
and one a 2;-dropout - at 2; ~ 6.5 they have redshifts in the transition range where small 
changes in flux in the bands can result in a change in their dropout designation. Given the 
similar redshifts, SEDs, and luminosities, and the quite small separation in projection, it 
is possible that the two sources are simply components of the same source with a physical 
separation of 2 kpc. We think it is unlikely they are multiple images of the same source since 
our lensing models do not predict such a split pair. Both sources show dual-core morphology 
(see §4.6), and their four components seem well aligned over a scale of 6" (see the J-band 
image in Fig. 4). We have kept them as separate entities throughout our discussion but 
suspect that it is highly likely that they are spatially associated. Given the small separation 
and similarities these objects will likely be components of a group that will ultimately merge. 



4.4. Surface Density of Bright High-Redshift Candidates 



We compare our samp le (Fig. 8) with the d istribution of known galaxy candidates at 
2; ~ 6, based on the LF of iBouwens et al.l (120061 ) and scaled to a search area of 25 square 
arcminutes. To a limiting magnitude of 25.6, a survey in an unlensed fleld is anticipated to 
flnd ~1.4 candidates at 2 ~ 6. We, however, flnd 4-5 such bright 2; ~ 6 — 7 candidates behind 
lensing clusters. This effectively illustrates how lensing from foreground clusters can help us 
identify high-redshift galaxies at relatively bright magnitudes - which is helpful for numerous 
follow-up studies, including spectroscopy. The observed enhancement is expected because 
of t he steep slope of the LF expected at higher luminosities {i.e., —d{\ogd(^)/d\ogL > 1: 



e.g 



Broadhurst et al.l (Il995l )). Such steep slopes are seen brighter than L* where there is an 
exponential cut-off, but not fainter than L*. For Al703-iDl, its intrinsic 2;-band magnitude 
is approximately 25 after correcting for the 3x model magnification factor (see Table 4). 
According to the 2; ~ 6 LF, one would expect to require 40-250 square arcminutes to find 
one such bright candidate. Since the discovery of Al703-iDl was made over an area of 
approximately 7 square arcminutes (after correcting for the 3x dilution factor), it would 
seem to be a rather serendipitous discovery. For the four other dropout candidates listed in 
Table 4, if we assume a nominal magnification factor of 9 in the cluster field, their intrinsic 
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magnitudes are AB 27 — 28. Over two square arcminutes in the source plane, the z ~ 6 
LF predicts approximately 2-5 such galaxies. The four such galaxies found are therefore 
consistent, within the uncertainties, with the known distribution of high-redshift galaxies. 



4.5. Rest-frame Optical Fluxes of Lower-Luminosity Galaxy Population 



Perhaps the most interesting property to examine in our small sample of highly mag- 
nified galaxy candidates is their stellar mass. The stellar masses tell us about the past star 
formation rate in a galaxy and therefore its likely contribution to the reionization of the 
universe as well as the production of metals. The stellar masses estimated in the previous 
section are useful for assessing the contribution of individual galaxies to the reionization 
of the universe. However, what we really want to know about the galaxy population at 
2; ~ 6 — 7 is their stellar mass density - since we can use this to constrain the SFR history at 
z > 1 . It is difficult to determine this directly from the above measurements since we only 
have five high-redshift galaxy candidates from which to derive a mass function (and also the 
mas ses we derive will be s ensitive to uncertainties in the model magnification factors: see, 
e.g. 



Bouwens et al.ll2009fcr ) 



A better way of estimating the total stellar mass density is to start with determinations 
of the UV LF at 2; ~ 6 — 7 and then use the measurements of UV-to-optical colors of 
2; ~ 6 — 7 galaxy candidates to estimate the stellar mass function of the galaxies (given the 
link between rest-frame optical colors and our estimates of stellar mass). Such colors have 
been available for the more UV lu minous galaxy candidates at 2 ~ 6 — 7 (lEyles et al.l 12005 



Yan et al.ll2006l : IStark et al.l 120071 ). but until the present have been lacking for the less UV 
luminous galaxy candidates at 2; ~ 6 — 7. This is where our sample of highly magnified 
galaxy candidates at z ~ 6 — 7 is useful since they provide a luminosity-to-mass calibration 
for the less luminous galaxies that have previously lacked colors. 

For the sake of uniformity, we will use the H — 3.6/im colors of our sample of high- 
redshift galaxy candidates to represent the rest-frame UV-to-optical colors. We plot this 
color distribution for the four strongly lensed galaxy candidates in our sample with IRAC 
fluxes in Fig. [91 We compare this distribution with the colors o f more luminous g alaxy 
candidates at 2 ~ 6 - 7 identified in the HUDF jYan et al.ll2005l : Ihabbe et al."1l2006L No 
significant differences are observed between the H — 3 .6/im colors of our sample and the 



more luminous galaxy candidates iLabbe et al. I (120061 ) examined in the HUDF. A t-test 



yields a significance of 51%, suggesting no large differences between the ages (and mass-to- 
light ratios) of the lower-luminosity galaxy population at high redshift and that at higher 
luminosities. This has obvious implications for the total stellar mass density at 2; ~ 6 — 7 
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and therefore the total star formation rate density at z > 7 (jChary II2008I ). since it allows 
us to estimate the fraction of the total stellar mass density at z ~ 6 — 7 that come from 
high- and low-luminosity systems. Our result does not suggest an increase in the number of 
low-luminosity galaxies. It is also consistent with the hypothesis that more massive galaxies 
are unresolved aggregates of lower-luminosity galaxies. 



4.6. Morphology of High-Redshift GalcLxies 

Strong lensing enables us to resolve otherwise unobservable spatial structures in high- 
redshift galaxies, and thereby investigate the structure of young galaxies and the physical 
processes that drive both their growth and vigorous star formation. The observed and in- 
ferred physical sizes and separations of the resolved components of the three sources reported 
here are given in Table |H For Al703-iDl, there are two knots separated by 0.45". Each 
of them can be fitted with a two-dimensional Gaussian with FWHM of 0.66" and 0.51", 
respectively. After demagnification (Table 4), this pair of components is separated by just 
0.9 kpc, suggesting a merging event that drives star formation. Both the candidates in 
CL0024-I-16, which themselves are probably related as discussed in §4.3, seem to consist of 
two components, which are best seen in the Jno image. 

It is interesting to note that the 2;-dropout in Abell 1689 (Al689-zDl) also has double 
nuclei, separated by ~ 1 kpc (Table 4), in the source plane. Double nuclei seem to be 



almost ubiquitous in this sample (see Table 4), and even A2218-iDl (IKneib et al.l l2004j ) 
appears to show significant substructure. One interpretation is that these galaxies may 
consist of multiple star-forming regions that are UV luminous. The separation between the 
two sub-structures seems to be sufficiently large (Tab le 4) that they r nay b e two merging 



systems. For A1689-zDl or A1703-iDl the models of iLimousin et all fl2007l . 120081 ) do not 
predict any counter images for these sources, and there is a lack of any luminous foreground 
sources that would be responsible for this lensing. Similarly as noted in §3, our models 
for CL0024-I-16 do not suggest multiple images for the sources. The complex and multiple 
nature of all these sources certainly suggests, as we would expect from numerical models, 
that interaction, merging and clumpy star formation are all aspects of the very dynamical 
nature of the galaxy build-up process at these redshifts. 
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SUMMARY 



We have found three bright, s trongly lensed i- an d 2;-dropouts in the fields of the 
massive clusters Abell 1703 (see also lRichard et al.ll2009l ) and CL0024+16. Their large color 
decrements make the selection of these high-redshift galaxy candidates robust (with the 
exception of CL0024-zDl for which the evidence is little weaker), and detailed cluster mass 
models enable us to derive the intrinsic properties of the high-redshift galaxies. The similar 
redshifts, SEDs, and luminosities, and the quite small separation in projection, of the two 
sources in CL0024+16 suggest that they may be spatially associated with a separation of 
~ 2 kpc. Using stellar-population models to fit the rest-frame UV and optical fiuxes, we 
derive best-fit stellar masses on the order of 10^~^° Mq. Even after correcting for the lensing 
magnification (a factor of 3), Al703-iDl remains intrinsically luminous. The three galaxy 
candidates are extended, with intrinsic size s between 0.2 and 2 kpc, and resolved into a 
pair of nuclei. They, along with Al689-zDl (IBradley et al.ll2008l ). which is also double with 
component separation of ~ 1 kpc in the source plane, suggest that merger driven growth 
and star formation may be a defining characteristic of high-redshift galaxies. Bright resolved 
galaxies such as these will be important targets for spectroscopy with the largest ground 
based telescopes and future large-aperture space telescopes like JWST. 
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Figure Captions 



Fig. 1. — Composite color (r625, ^775, ;Z85o) ACS image of Abell 1703. The z = 6 critical curves 
(magnification factor fi > 75) for the lensing model of lLimousin et al.l (120081 ) are overplotted. 
The positions of the shallow (one orbit each) NICM0S/NIC3 FllOW and F160W images 
are illustrated with the four blue boxes and an orange box, respectively. The location of 
Al703-iDl is marked with a red circle. 

Fig. 2. — Composite color (r625, ^775, -2350) ACS image of CL0024+16. Th e z = 6.5 critical 
curves (magnification fi > 200) for the lensing model of IZitrin et al.l (120091 ) are overplotted. 
The blue boxes show the regions wit h deep NICM0S/NIC 3 FllOW imaging observations. 
Two of these NIC3 pointings are from lRichard et al.l (120081 ) and two are from HST program 
GO-10874. The positions of the two dropout galaxies CL0024-iDl and CL0024-zDl are 
marked with the red circles. 



Fig. 3. — Cutout images of the z-dropout galaxy Al703-iDl. Each cutout is 15"on a side. 
The source is detected in the ACS z band and NICMOS FllOW and F160W bands, but not 
in the five optical bands at shorter wavelengths. Th e image formed from the F435W, 
F475W, F555W, F625W, and F775W optical bands jSzalav et al.lll999h is presented here 
and labeled "Det." The J-band data taken with the Gemini North telescope are also shown. 
This source is clearly resolved into two extended nuclei. With an iJ-band magnitude of 23.9, 
it is one of the most luminous galaxy candidates at 2; ~ 6 ever found. 

Fig. 4. — Cutout images of the i-dropout galaxies CL0024-iDl (marked with a red circle) 
and CL0024-zDl (green circle). Each cutout is 15"on a side. T he image forme d from 
the F435W, F475W, F555W, F625W, and F775W optical bands JSzalav et al. l ll999h is also 
presented here and labeled "Det." CL0024-iDl is detected in the ACS z-band, NICMOS 
FllOW and F160W bands, and the IRAC 3.6/im and 4.5/im bands while CL0024-zDl is 
detected in the NICMOS F110W,F160W and the IRAC 3.6/im band. Neither source is 
detected in the five ACS optical bands at shorter wavelengths. The J-band data taken with 
the Walter Baade (Magellan) telescope are also shown. In the Jno band, each source displays 
a dual-core structure, and the four components are well aligned. 

Fig. 5. — Best-fit stellar-population models for A1703-iDl. The stellar masses are intrin- 
sic values, corrected f or the magnification factor of 3 as predicted by the lensing model of 
Limousin et al.l (120081 ). The vertical bars denote the la flux uncertainties and the horizontal 
bars represent the filter bandwidths. Non-detection in a given band is marked with a 2a 
upper limit in flux and a down arrow. The deep ACS data eliminate lower-redshift sources. 
The red and blue curves represent the z > 5 fits with solar metallicity (Z=0.02) and sub- 



- 19 - 



solar (Z0/5O) metallicity, respectively, and the cyan curve for the z = 5.827 fits with solar 
metallicity. The green and brown curves are for z < A fits and have high values of xt- 

Fig. 6. — Best-fit stellar-population models for CL0024-iDl. The stellar masses are intrinsic 
values, corr e cted f or the magnification factor of 6.3 as predicted by the lensing model of 



Zitrin et al.l (120091 ). The green and brown curves are for z < 4 fits and have high values of 



xl- See the caption of Fig [51 

Fig. 7. — Best-fit stellar-population models for CL0024-zDl. The stellar masses are correcte d 
for the magnification factors of 6.2 as predicted by the lensing model of IZitrin et al.l (120091 ). 
The z > 5 fit gives a much lower value oi xl ^ ^ than for the z < 4 fit where xl ~ 2.5. Also 
see the captions of Fig [5] and Fig [61 

Fig. 8. — i-dropout surface density. The hatched are a shows the nurnber o f z-dropouts per 
20 square arcminutes in the field, as determined by iBouwens et al.l (120071 ). The observed 
brightness for the five dropout galaxies in Table 2 are marked with stars, while their intrinsic 
(demagnified) if-band magnitudes are marked with the solid circles. The high magnification 
results in a number of bright sources (though they are intrinsically low luminosity, in general). 
As suggested in a number of studies, cluster magnification can result in bright z ~ 6 galaxies 
{H < 25) that are rare in the field (see §4.3). 

Fig. 9. — The H — 3.6/xm color distribution (red histogram) for the lower luminosity, highly 
magnified galaxy candidates at z ~ 6.5 — 7.6 in our selection (Table [2]). This distribution 
should provide an approximate measure of the size of the Balmer break and therefore provide 
us with an approximate measure of its age (though this color will also depend upon the dust 
content). In general, redder galaxies would be expected to be older and bluer galaxies 
younger. Also shown is the color distribution for the more lu r ainous galaxy cand idates 
at 2 ~ 6.7 - 7.2 in the HUDF (black histogram: IVan et al.l[2005l : Ihabbe et aini2006h . Our 
sample of lensed galaxies has similar H — 3.6/im colors to those of the more luminous sample. 
This suggests that the lower-luminosity galaxy population may have similar mass-to-light 
ratios as higher-luminosity galaxies (see §4.4). 
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Table 1. HST NICMOS Data around Massive Galaxy Clusters Used to Search for i— and 

2;— Dropout Galaxies at z ~ 6 — 7 



Area 5cr Depth AB Magnitude^ 
Name (arcmin^) 2775*^ Zgso Jno -f^ieo ReF 

New i- and z-Dropout Search Data (HST GO 10874) 



CL0024+16 


1.5 


27.8 


27.4 


26.7 


26 


6 


[1] 


CL1056-03 


2.2 


27.5 


27.3 


26.6 


26 


5 


[1] 




i-Dropout Search Dat 


d 








MS1358+61 


1.4 


28.0 


27.5 


26.8 


26 


7 


[2] 


CL0024+16 


1.4 


27.8 


27.4 


26.8 


26 


7 


[2] 


Abell 2218 


1.2 


27.9 


27.6 


26.8 


26 


7 


[2] 


Abell 2219 


1.4 


27.5 


27.2 


26.8 


26 


7 


[2] 


Abell 2390 


1.5 


26.6 


26.9 


26.8 


26 


7 


[2] 


Abell 2667 


1.5 


26.5 


26.9 


26.8 


26 


7 


[2] 


Abell 1689 


5.7 


28.1 


26.8 


26.4 






[3] 


Abell 1703 


3.1 


27.9 


26.7 


26.0 






[4] 


lE0657-56^ 


2.7 


27.6 


27.2 


26.0 






[4] 


Abell 1835 


0.7 


27.1 


27.2 




26 


7 


[4] 


AC114 


0.7 


27.1 


27.2 




26 


7 


[4] 



^5cr limits assume a 0.3"-diameter aperture for ACS/WFPC2 and 0.6"-diameter aperture 
for NICMOS. The source photometry in this paper uses a larger aperture, therefore leading 
to shallower limiting magnitudes. 

^The depth of the deepest single-band optical (< 0.8/im) image available over the cluster. 
Note that some of the clusters listed here have very deep data in two or more bands (par- 
ticularly for MS1358+61, CL0024+16, Abell 2218, Abell 1689, Abell 1703), so the effective 
depth of the combined optical data is often >0.4 magnitude deeper than tabulated here. 



Bradlev et all (120081 ). [4] iBouwens et all (l2009a[ ). 



References : [1] This paper. [2] iRichard et all fl2008h . [3] iBroadhurst et all (120051 ): 



^These data were previously used by lBouwens et al.l (l2009al ) to search for z and J dropout 
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galaxies. We use these data to search for bright i dropouts. 
^The "Bullet" cluster. 



Table 2. Photometry of Bright z > 6 Dropout Objects in Cluster Fields (AB Magnitude) 



Name aj2ooo ^J2ooo ^^625 «775 ^sso -^iio -^^leo 3.6/im 4.5/im 



Al703-iDP 


13 


15 01.41 


+51 


48 


25.9 


> 


28.2 


26.J 


>±0.3 


24.2 


± 


0.1 


24.0 


± 


0.1 


23.9 


± 


0.1 


23.1 


± 


0.3 


23.5 


± 


0.4 


CL0024-iDl^ 


00 


26 37.78 


+17 


10 


40.0 


> 


28.1 


> 


27.9 


26.0 


± 


0.2 


25.1 


± 


0.1 


25.0 


± 


0.1 


24.4 


± 


0.2 


24.4 


± 


0.3 


CL0024-zDl^ 


00 


26 37.93 


+ 17 


10 


39.0 


> 


28.1 


> 


27.9 


27.3 


± 


0.8 


26.0 


± 


0.2 


25.6 


± 


0.2 


24.5 


± 


0.5 


24.8 


± 


0.5 


Al689-zDl^ 


13 


11 29.96 


-01 


19 


18.7 


> 


27.8 


> 


27.8 


> 


27 


5 


25.3 


± 


0.1 


24.7 


± 


0.1 


24.2 


± 


0.3 


23.9 


± 


0.3 


A2218-iDr''^ 


16 


35 54.40 


+66 


12 


32.8 


> 


27.6 


> 


27.2 


25.1 


± 


0.1 


24.3 


± 


0.1 


24.1 


± 


0.1 


23.7 


± 


0.3 


23.9 


± 


0.3 



^This paper. These two sources are only 2.5" apart, have similar redshifts and are probably spatially associated (See Table 
3 and 4 and the text). 



Bradley et al.l ((20081). 



00 



Egami et al.l (120051 ). The r- and i-band magnitudes are from our own measurements. 



'^Component b. The other component (a) has a similar magnitude, but its IRAC fluxes are unavailable because of contam- 
ination from a neighboring source. 
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Table 3. Best-fit Stellar- Population Model Results 



Object 


^phot ^form 


Mass*' Age'' 
(109 Mo) (Myr) 


Av 






High Redshift, Solar Metallicity {Z 


= 0.02) 




Al703-iDl 

Al703-iDl 
CL0024-iDl 
CL0024-ZD1 


5.9 ±0.1 6.2 
5.827 6.2 

6.4 ±0.1 6.7 

6.5 ±0.3 7.1 


9.2 ± 3.2 55^^^ 
12 ± 0.6 64^^° 
1.5 ±0.1 43^}^ 
2.0 ± 0.4 811^^ 


0.0 ± 0.2 
0.0 ±0.1 

0.0 ±0.1 
0.0 ±0.1 


1.2 
1.1 
0.7 
0.1 


High Redshift, Subsolar Metallicity (Z 


= 0.0004) 




Al703-iDl 
Al703-iDl 
CL0024-iDl 
CL0024-ZD1 


6.0 ±0.1 6.5 
5.827 5.9 
6.5 ±0.1 7.0 
6.7 ±0.4 7.4 


12 ± 4.6 64+^^ 

13 ± 3.4 13tl% 
2.1 ±0.8 72;^4 
2.5 ± 0.2 9Qtl7^ 


0.2 ±0.2 
1.2 ±0.1 
0.0 ±0.2 
0.2 ±0.4 


1.1 
1.7 
0.5 
0.1 


Low Redshift, Solar Metallicity (Z -- 


= 0.02) 




Al703-iDl 
CL0024-iDl 
CL0024-ZD1 


1.2 ±0.1 1.4 

1.3 ±0.2 1.5 

1.4 ±0.5 1.5 


2.4 ± 0.4 4541^^^ 
0.4 ±0.1 571^^1 
0.4 ± 0.2 5711^^5 


0.8 ±0.3 
0.3 ±0.4 
0.7 ±0.9 


9.3 
2.8 
0.6 



Low Redshift, Subsolar Metallicity (Z = 0.0004) 



Al703-iDl 1.2 ±0.2 1.4 4.2 ± 0.9 454+^f 1.8 ± 0.3 11.4 
CL0024-iDl 1.3 ±0.2 1.5 0.7 ± 0.2 1.4 ± 0.4 4.2 

CL0024-ZD1 1.4 ±0.6 1.6 0.6 ± 0.4 4041^°^ ig±iQ o.9 
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^Formation redshift calculated from the fitted redshift and SSP stellar 
age. 

^Best-fit stellar mass, corrected for the cluster magnification. 
^Age of the single- burst stellar population (SSP model). 



Table 4. Intrinsic Properties of Bright Dropout Objects 



Name 


Photometric 


Observed Size Demagnification 


Luminosity 


Size 


Component Relative 




Redshift 


(") Factor 


(10'^^ ergs s"-^) 


(kpc) 


Separation (kpc) Flux 



Al703-iDl 

CL0024-iDl 

CL0024-ZD1 

A1689-ZD1 

A2218-iDl^ 



5.95 ±0.15 
6.45 ±0.15 
6.55 ±0.35 
7.6 ±0.4 
6.65 ±0.1 



0.66 ± 0.5P 
1.1 ± 0.53^^ 
0.6 ± 0.5^ 
1.6 ± 1.0^ 
3.0 X 0.3^ 



3 

6.3 
6.2 
9.3 
15 X 1.7^ 



21 
5.5 
4.8 
6.9 
3.2 



2.2 ± 1.7^ 
1.0 ± 0.5^ 
0.5 ± 0.4^ 
0.9 ± 0.5^ 
0.2 



0.9 
1.3 
1.4 
1.6 



1.0/0.57 
1.0/0.09^^ 
1.0/0.51^ 
1.0/0.62^ 



Tor two separate nuclei. 



Egami et al.l (120051 ). These two sources are only 2.5" apart, have similar redshifts and are probably spatially associated (^ee 



Table 2 and 3 and the text). 

^Along the major and minor axes, as the image is greatly stretched due to a gravitational lensing effect. 
"^Measured in the J band. 
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